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Abstract 

Phase space representations of the dynamics of the quantal and classical cat 
map are used to explore quantum-classical correspondence in a K-system: as 
7i — i- 0, the classical chaotic behavior is shown to emerge smoothly and exactly. 
The quantum dynamics near the classical limit displays both exponential sep- 
aration of adjacent distributions and long time relaxation, two characteristic 

features of classical chaotic motion. 
PACS numbers: 05.45, 03.65.S, 05.40 
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Understanding the correspondence principle for bound conservative systems which are 
classically chaotic presents a serious challenge |Ij]. The classical system is characterized by 
positive K-entropy which is reflected, for example, in the exponential divergence of initially 
nearby trajectories and in the relaxation || of distributions in the long-term limit. By 
contrast, the quantum dynamics is not even ergodic, and displays long time recurrences. 
Current views on correspondence in conservative chaotic systems range from simple belief in 
the correspondence principle, to reliance on decoherence induced by coupling to an environ- 
ment H, to calls for the overthrow of quantum mechanics due to its purported inability to 
display chaos in the correspondence limit Resolving this controversy requires studies of 
quantum dynamics close to the classical limit, a task thus far made impossible by numerical 
difficulties as % — > 0. 

In this Letter we demonstrate, using numerical results based on an exact quantum prop- 
agator 0] , the approach to the classical limit of the chaotic Arnold cat || . The nonsingular 
character of our phase space propagator as h — > allows direct insight into correspondence 
in this bounded chaotic system. In particular, we show that quantum mechanics does dis- 
play chaos as % — > 0, including exponential divergence of initially nearby distributions and 
long time relaxation. The agreement between quantum and classical dynamics extends over 
longer times as % decreases ||. 

Classical-quantum correspondence is best examined |7|] in the Liouville picture classically 
and in some phase space representation of the density matrix (e.g. the Wigner representation 
H) quantum mechanically. This approach provides an overlap of objects and concepts which 
does not exist if one attempts to utilize, for example, classical trajectories and quantum 
wavef unctions. The advantages of this approach are evident in the results below where 
operationally chaotic behaviour is demonstrated in quantum mechanics. 

The classical cat map is a K-system which corresponds to the classical dynamics of a 
kicked oscillator with Hamiltonian || 

oo 

H = p 2 /2fi + ex 2 /2 8(s-t/T). (1) 

s=— oo 
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restricted to a torus < x < a, < p < b. The proper quantization of this system and its 
approach to the classical limit has been extensively considered PJIU|1| . Early quantizations 



|TjJ of this system were unsettling insofar as: (a) only a restricted class of classical tori (those 
with ab = hN, N integer) could be quantized; (b) the solutions violated the uncertainty 
principle and (c) the quantum propagator thus obtained did not reduce to the classically 
chaotic cat map in the H — > limit. 

A recently introduced @] phase space quantization eliminated these problems and yielded 
an analytic quantum propagator. Specifically, the one time-step quantum propagation of 
the Wigner-Weyl representation of the quantal density operator p®(x,p;t) in phase space 
(x, p) is given by 

1 oo 

pQ( X ,p;T) = - J2 Plm(T)fn, m (2) 
0,0 n,rrt=-oo 

where 

oo 

pl m (T)= £ G(n,m;k,l)pl(0). (3) 

k,l=— oo 

The propagator G of the Fourier coefficients is given by 

G(n, m; k, I) = e ina{kl - nm) du ^ du' e- i7r ^ u+an)mod nV^^/a . 

Jo Jo 

^2ni(l— ra)v ^— i7r-^[(i^'+a/)mod l] 2 /a ^iirrjv' 2 /a ^2iri(k—n)v' /^"j 

where r] = Tb/pa G Z (the set of integers) and £ = —eTa/b G Z obtain in scaling the 
original Hamiltonian to the unit torus, a = h/ab acts as a dimensionless form of Planck's 
constant for this problem and f n ,miPi x ) — exp{27ri (rzp/a + mx/b)} is the Fourier expansion 
basis. The Arnold cat map corresponds to the choice rj = £ = 1. Equation (Q) can be 
integrated exactly [|ll[] to provide an analytic propagator for all H. This result contrasts with 
the classical propagator given by Eq.(^) but with p® m (T) replaced by the classical form 

oo 

p n ,m{T)= Gc{n,m;k,l)p k> i(0) 

k,l=— oo 

G c (n, m\ k, I) = S^jr.fam). (5) 



Here T • (n, m) is the vector resulting from premultiplying the the vector (n, m) with the 
transpose (f> T of the matrix governing the point dynamics: 
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with x = x/a,p = p/b. 

Note that, by contrast with the classical case, the quantum propagator G mixes contribu- 
tions from all k, I components to produce each n, m. It shares with the classical propagator 
a fixed point at (n, m) = (0,0). However, it also has an infinite number of additional fixed 
points for rational a = u/v; these lie in the quantal Fourier space at all (wv,zv) with 
u,v,w,z G Z. These fixed points affect the quantal evolution of a distribution in the fol- 
lowing ways: (a) near the classical limit, (a C 1 ort; > m > 1), these fixed points are far 
out in the Fourier space, and hence affect structure at extremely small scales. Thus, it is in 
the highly oscillatory regime (whether the structure is prepared initially or acquired during 
the evolution) that the dynamics will exhibit quantal evolution distinct from the classical 
dynamics; (b) at larger a, these fixed points are both more dense in the Fourier space and 
exist at smaller values of (n,m), indicating a far more substantial deviation from the classi- 
cal evolution; (c) even though there are no such demonstrable fixed points for irrational a, 
the propagator is a smooth function of a; the qualitative conclusions from these arguments 
should therefore extend to all a. 

The evolution of distributions which are initially of the Gaussian form 

N x exp[-(x - \x ) 2 /\ 2 a 2 x ] exp[-(p - Ap ) 2 /AVJ] (7) 

is displayed below. Here N\ is a normalization factor, (\x , \p ) specifies the location of 
the Gaussian of width (Xa x , \<r p ) in a phase space of dimension \ 2 ab. The variable scaling 
factor A can be used to approach the classical limit. That is, as A increases the volume of 
phase space \ 2 ab increases and a = h/\ 2 ab, the effective Planck's constant of the problem, 
decreases. Hence, by increasing A we approach the classical limit while preserving the ratio 
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(\ 2 a x <T p / X 2 ab) of the volume of the initial distribution to the volume of phase space. We 
note that we work in units such that h — 1. 



Figures 1-3 show gray-scale contour maps of the evolution [|12[] of the Wigner function 
for a = 10 _1 , 10~ 2 and 10 -5 . To simplify the presentation we report the absolute value of 
the Wigner representation of the density. The Wigner-Weyl function for a = 10~ 5 is always 
positive definite and numerically identical to the classical Liouville density over the indicated 
time scale. We note that the classical dynamics is scale-invariant; hence the classical phase- 
space figures are independent of A and are the same for all three cases and identical to the 
a = 10~ 5 case. The analytic results |4| ensure that the a — > limit is indeed classical 
mechanics but these computations allow a direct examination of the approach to this limit. 

For a = 0.1 (Fig. 1) the system is seen to be in the deep quantum regime. The propa- 
gating wavepacket shows only the barest resemblance to the classical result and interference 
structures dominate. There is no relaxation of the dynamics over the indicated timescale 
and the expectation value of operators (not shown) oscillate with large amplitude. Although 
not evident in this figure, the Wigner distribution has many positive and negative regions. 
For a = 0.01 (Fig. 2) the classical density striations begin to appear, although still accom- 
panied by quantum effects. Various parts of the distributions appear to interfere with one 
another giving interference fringes. Once again relaxation is not evident after eight time 
steps, although a large coarse graining will result in a uniform phase space distribution. 
Separate calculations show that the expectation value of operators in this a regime display 
small fluctuations about the equilibrium value. Finally, as noted above, the a = 10~ 5 result 
agrees exactly with the figures from the classical evolution: The distribution has relaxed || 
even over the short time-scale shown. The approach to classical mechanics is clearly smooth 
and devoid of singularities. 

Relaxation is but one characteristic of chaotic motion, evident at long times. Classical 
chaos, however, is typically manifested in the short time instability associated with exponen- 
tial divergence of adjacent trajectories. Standard arguments (T3| imply that such exponential 
separation is not possible in quantum mechanics. Here we demonstrate that this is not the 



case; we have obtained similar exponential divergence for compact phase space distributions, 
both classically and quantally close to the classical limit. In the absence of point trajectories 
and an appropriate tangent space, we examine the time dependence of the distance A be- 
tween two nearby, initially well localized Gaussians |TJ|]; the results are shown in Fig. 4. Also 



shown for comparison is a straight line with slope dictated by the K-entropy of the classical 
cat map; since the system is uniformly hyperbolic, this is exactly equal to the Lyapunov 
exponent of the underlying point dynamics. The distance in phase space between the two 
distributions is determined by the centroids: A = [({x) 1 — (x) 2 ) 2 + ({p) 1 — (p) 2 ) 2 } 1 ^ 2 , where 
() f denotes the average over the i th distribution (i = 1, 2). The quantity A(i) was evaluated 
both classically and quantum mechanically for various values of a, b. 

The time dependence of A was found to be essentially identical in classically and quan- 
tally propagated distributions over times up to eight steps, for the a, b (and hence a) in 
the range shown. At early times, dependent on a, A increases essentially exponentially [see 
Fig. 4], with increasing adherence to the classical K-entropy line as a is reduced. At times 
longer than that shown in Fig. 4 A tends to decrease in a generally similar way quantum 
mechanically and classically, a consequence of saturation [JHjJ and relaxation in the finite 
size phase space. When a is increased, since there is no relaxation, A continues to oscillate 
at larger times quantum mechanically, whereas the classical dynamics relaxes. The com- 
bination of short-time exponential divergence and long-time relaxation demonstrates the 
loss of information characteristic of chaos for both the classical and near-classical quantal 
dynamics. 

Finally, we note that the availability of an analytic quantum propagator for a K-system 
allows a study of the expected time-dependence of classical-quantum agreement as a func- 
tion of h. Preliminary results suggest that increasing the value of a does, in general, decrease 
the time-scale over which quantum and classical dynamics agree. However, the extent of 
agreement also depends strongly on various other factors, including the nature of the ini- 
tial distribution. The variation within the quantum-classical deviations for different initial 
conditions at a given value of a and t spans many orders of magnitude. Further studies 



on this issue await our application of our general quantization approach to classical maps 
which have smaller K-entropy than the classical cat map. That is, the large value of the K- 
entropy for the cat map forces relaxation of the dynamics over a rather short time. Further, 
it also rapidly saturates the exponential divergence [Fig. 4] which could otherwise extend to 
significantly longer times. 

In summary, although the literature suggests that demonstrating quantum-classical cor- 
respondence becomes increasingly difficult (if not impossible JyJI) with increasing classical 
chaos, we have been able to show this correspondence in a classical K-system. Secondly, 
we have shown how correspondence emerges smoothly from the quantal behavior as a — > 0. 
These results make clear that the dynamics of distributions in the non-ergodic quantal cat 
map display the essential attributes of the classical K-system as the system approaches the 
classical limit. They constitute the first demonstration of the emergence of characteristi- 
cally classically chaotic behavior from quantum mechanically non-ergodic behavior as one 
approaches this limit. 
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FIGURES 

FIG. 1. Snapshots of \p®(x,p, t)\. Densities are represented by five shadings of grey denoting 
five evenly spaced probabilites between 7 x 1CT 5 to 1 x 10~ 6 . 

(a) t = 0,(b) t = 2T, (c) t = 4T,(d) t = 8T; a = 0.1 (Xa = Xb = 3.162). 

FIG. 2. As in Figure 1 but with a = 1(T 2 {Xa = Xb = 10). 

FIG. 3. As in Figure 1 but with a = 1(T 5 (Xa = Xb = 316.2). 

FIG. 4. Short-time exponential separation of centroids for a = b = 50.0 (a = 4x 10 -4 ) (dashed) 
; a = b = 28.79 (a = 1.207 x 10~ 3 ) (solid) ; a = b = 16.57 (a = 3.64 x 10~ 2 ) (dotted). See text for 
other parameters. Similar results are obtained by propagating the same distributions classically. 
The straight line with a slope of 0.9624 shows classical exponential separation of adjacent point 
trajectories. 



10 



